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We examine the robustness of single-molecule tunneling diodes to thermal-environmental effects. The diode 
comprises three fragments: two different conjugated chemical groups at the boundaries, and a saturated moiety 
in between, breaking conjugation. In this setup, molecular electronic levels localized on the conjugated groups 
independently shift with applied bias. While in the forward polarity a resonance condition is met, enhancing 
conductance, in the reversed direction molecular electronic states shift away from each other, resulting in 
small tunneling currents. In the absence of interactions with a thermal environment (consisting e.g. internal 
vibrations, solvent), rectification ratios reach three orders of magnitude. We introduce decoherence and 
inelastic-dissipative effects phenomenologically, by using the “voltage probe” approach. We find that when 
Id ^ V, with 7d the interaction energy of electrons with the environment and v the tunneling energy across 
the saturated link, the diode is still highly effective, though rectification ratios are cut down by a factor of 
2-4 compared to the coherent limit. To further enhance rectification ratios in molecular diodes we suggest a 
refined design involving four orbitals, with a pair of closely spaced states at each conjugated moiety. 


I. INTRODUCTION 

One of the simplest building blocks in ordinary semi¬ 
conductor circuitry, essential for realizing molecular- 
based electronics, is the rectifier (diode). In fact, the 
Aviram-Ratner proposal for an organic molecular rec¬ 
tifier largely initiated the field of molecular electronics 
in 1974i. Since then, molecular rectifiers have been 
fabricated from thin films, self-assembled monolayers 
(SAMs)^“— , and single molecules^Ti^. However, in general 
reported rectification ratios were rather low. Only re¬ 
cently have SAMsiSr— and single molecule junctionsi^d^ 
provided robust molecular rectifiers with high conduc¬ 
tance, and rectification ratios of two and even three or¬ 
ders of magnitude. 

What is the mechanism of electron current rectifica¬ 
tion in molecular junctions? The literature includes sev¬ 
eral constructions based on molecules with an asymmet¬ 
ric backbone, different molecule-metal linkers at the two 
ends, or electrodes of different materialsi^. Fundamen¬ 
tally, from the principles of quantum transport, it is ap¬ 
parent that coherent conduction of noninteracting elec¬ 
trons, as described in scattering theory by the Landauer 
formula, cannot materialize rectification when energy lev¬ 
els are fixed, independent of applied biasii. Two neces¬ 
sary conditions should be simultaneously met for achiev¬ 
ing the diode effect: many body interactions should play 
a role, and a spatial asymmetry should be built into the 
junction, see e.g. Refsi^^d^, exemplifying charge and heat 
rectification. 

In many proposals, however, the role of many body 
effects is not spelled out, and it is only included at 
the level of mean-field, by assuming a screening inter¬ 
action which creates a voltage drop across the molecular 
bridgeiSi2£. In particular, theoretical and computational 
investigations of tunneling diodes employ molecules that 
are asymmetrically coupled to the electrodes. Using 
model calculations or first-principle simulations, it is as¬ 
sumed or demonstrated that the junction develops an 
asymmetric potential profile under bias in the forward 


and backward direction, see e.g. Refsi ^^’^^~ — . This in 
turn results in an asymmetric IV characteristics. The 
development of a voltage drop inside the molecule, how¬ 
ever, can be only justified when many body effects (es¬ 
sentially electron-electron interactions) play a role in the 

junctional. 


In this paper, we are concerned with single molecule 
rectifiers made of three fragments, e.g., donor-cr-acceptor 
junctions. While in their original proposal Aviram and 
Ratner had assumed vibrationally-assisted intramolec¬ 
ular electron transfer—, a related diode can be real¬ 
ized based on quantum tunneling, as presented in Fig¬ 
ure [T] and discussed in the next section. Tunneling 
diodes were investigated for example in Ref— 1, where 
a DFT -I- NEGF approach was employed for the study 
of coherent transport in a three-fragment l,2-bis(4- 
(phenylethynyl)phenyl), showing a rectification ratio of 
R ^ 1000. A similar approach (with a different back¬ 
bone) was adopted in Ref— 1, manifesting i? ~ 100. These 
studies and similar proposals employ the Landauer for¬ 
mula for coherent conduction, with the additional as¬ 
sumption of level shift under bias. 


The objective of the present work is to examine the 
robustness of tunneling diodes against environmental ef¬ 
fects, when electrons on the junction suffer from decoher¬ 
ence, inelastic effects, and energy dissipation due to their 
interaction with an “environment”, consisting of degrees 
of freedom such as molecular vibrations, the solvent, and 
other electrons. Naturally, since the rectification mecha¬ 
nism in tunneling diodes, depicted in Fig. [U relies on the 
sharp contrast between resonant transmission and deep 
tunneling, rectification ratios should deteriorate when in¬ 
coherent processes play a role. However, an estimate of 
this effect for realistic molecular diodes has not yet been 
provided. 
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FIG. 1. Illustration of the examined tunneling diode. Displayed are levels (a) at zero bias, (b) forward bias, providing a 
resonant situation beneficial for conductance, and (c) under reverse bias when levels are shifted away from each other resulting 
in small tunneling currents. 


II. MODEL AND METHOD 


We focus on recent proposals for highly rectifying tun¬ 
neling diodes built from three-fragment molecular back¬ 
bones such as l,2-bis(4-(phenylethynyl)phenyl)ethane 
junctiona^ij see Fig. [2j This asymmetric molecule in¬ 
cludes conjugated groups at each end of the junction and 
a weak tunneling link in the form of a saturated carbon 
bridge, breaking conjugation. IV characteristics of this 
system, as calculated with a DFT -|- NEGF approach, 
were excellently matched with the two-site (HOMO and 
HOMO-1) model Hamiltonian 


H = 


ei -h ^aeV v 

V €2 — \aeV ’ 


with £12 as the conjugated groups’ site energies at zero 
bias, V a small tunneling energy between the units, and a 
a phenomenological parameter which describes the frac¬ 
tion of the voltage drop within the molecule, the response 
to the screening interaction. At zero bias, ci = 0 and 
£2 < 0. Under forward bias defined as U > 0, a reso¬ 
nance situation can be met once £i -I- \aeV « £2 — \aeV^ 
with the shifted levels buried within the bias window, re¬ 
sulting in high currents. Under the reversed operation 
only tunneling currents contribute in coherent scenarios, 
see Fig. [T] Coherent transport can be described by the 
Landauer formalismi^. We construct the metal-molecule 
hybridization matrices with coupling strengths 7 l,_r, 


temperature T = eV = — p./;, e is the elec¬ 

tron charge. In Fig. [2] we display the IV characteristics 
and the corresponding rectification ratio 


R{y) = 


|/(^>0)l 

|/(U<0)| 


(4) 


for the l,2-bis(4-(phenylethynyl)phenyl)ethane molecule, 
reproducing rectihcation ratios up to i? ^ 1000 in the 
coherent limit, as discussed in Refi^. 

We now assess the role of a thermal environment on 
this predicted excellent rectifying behavior. The cou¬ 
pling of electrons to other degrees of freedom, e.g., 
internal molecular vibrations, solvent, other electrons, 
causes phase decoherence of transmitted electrons, in¬ 
elastic scattering processes, and energy relaxation. A 
detailed and accurate description of such processes in 
molecular conduction is a challenging task^, thus for the 
sake of simplicity and generality, we emulate the envi¬ 
ronment in a phenomenological manner using Biittiker’s 
probe approach^Si^. We mimic elastic and inelastic scat¬ 
tering of electrons on the molecule by coupling them to 
fictitious reservoirs, termed “probes”. Specifically, the 
two molecular sites are coupled to independent (ficti¬ 
tious) metals, identified by ’1’ and ’2’. These metals can 
exchange electrons with the junction, and their proper¬ 
ties are determined self-consistently to introduce relevant 
thermal effects. In direct analogy with Eq. m , we con¬ 
struct the molecule-probes hybridization matrices as 



■ 1 o' 


'0 o' 

7L 

0 0 

II 

0 1 


and organize the transmission function 

TlA^) = Tr[fLG'-(£)fflG“(£)], (2) 

defined in terms of the retarded and advanced Greens 
functions, G’’(£) = [el — H + ^{Tl + f^)]”^, G“(£) = 
[G'-(£)]t, I is the identity matrix. The current leaving 
the L terminal is calculated from 
2e /'°° 

II = -^ J c1£7l,_r(£)[/l(£) - /fl(e)], (3) 



'10' 


'0 0' 

Id 

0 0 

, r2 = 7d 

0 1 


with 7 d the electron-environment coupling strength. The 
essential convenience of the probe method lies in its di¬ 
rect correspondence to coherent transport: The current 
crossing the junction is calculated using the Landauer- 
Biittiker formula, by generalizing Eq. m to the multi¬ 
terminal case, 

O poo 

^^ = -h T. / derL,,(£)[/i(£)-/,(£)], (6) 


where /[/(c) = -|- I] ^ is the Fermi function at with7L,j(£) = Tr[f iG’'(£)f jG“(£)] and the Green’s func- 

the V = L,R terminal with the chemical potential and tions generalized to include the molecule-probe coupling. 
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G;(e) = \ei-H+ + fi + r^)]-^ G“(e) = 

[G'’(e)]^. Decoherence and energy exchange processes 
are implemented with “voltage probes”—, by demanding 
zero net charge flow into each probe, 



h.2 = 0 , 


(7) 


with e.g. the probe -1 current 


h 


Oo poo 

-r E / *'ri.,(e)[/i(e)-/,(e)]. 


( 8 ) 


The functions /i, 2 (e) assume a Fermi-function form at 
temperature T (uniform across the junction), and we 
search for ^ 1^2 solving Eq. O- Far-from-equilibrium, 
when the applied bias is large, an exact analytic solu¬ 
tion for this problem is generally missing, though unique¬ 
ness is guaranteed^^. We thus retract to a fully nu¬ 
merical procedure, recently implemented for electronic 
conductions^ and phononic heat transpor t^^i^^ , and em¬ 
ploy the Newton-Raphson method by iterating according 
to 





■ dh 

dh' 

-1 

T 

in 

dl2 


h 

h 

_ dfj,i 


k 



(9) 


converging to the unique set of roots for the probes’ 
chemical potentials. The voltage probe condition ([7]) al¬ 
lows phase loss processes and energy exchange. Decoher¬ 
ence without relaxation can be implemented as well with 
“dephasing probes”, by demanding the energy resolved 
currents at each probe to vanish, /i, 2 (e) = 0. These 
strict constraints translate into a set of linear equations 
for /i, 2 (e), which can be readily solved. 

The probe approach for decoherence and relaxation 
has found many applications in mesoscopic physics for 
understanding the phenomenology of charge^ 2 i 2 Si 21 and 
hea t^^i^^i^^i^^ transfer. More recently, several studies had 
implemented probes in molecular electronic problems, al¬ 
beit limited to linear response situations'^— or to the 
simpler dephasing probe approach^. 

Below we examine the functionality of tunneling diodes 
under inelastic effects by implementing the voltage probe 
approach. Setting the applied voltage, we iterate Eq. (0 
to convergence, then calculate the current in the junction 
from Eq. ([S]). Convergence is determined by ensuring 
that the probe potentials /ri _2 are evolving monotonically, 
IMiV ~ Ml, 2 1 ^ IMi ,2 ~ and by confirming that 

the leakage current to each probe is small, |/i, 2 / 7 l| < 
10“®. Our procedure generally required tens to hundreds 
of cycles to converge, with initial conditions at each bias 
prepared from converged results at lower voltages. 



FIG. 2. Tunneling diodes under inelastic ef¬ 

fects with parameters corresponding to the l,2-bis(4- 
(phenylethynyl)phenyl)ethane molecule (top panel), serving 
as the backbone of the diode, (a) IV characteristics of the 
junction, (b) rectification ratio R, and (c) level diagram show¬ 
ing the chemical potentials Rl,r at the two leads (dashed 
lines) and the eigenenergies of H, Eq. ([T]) (full lines). Param¬ 
eters here and throughout the paper correspond to a two-state 
model of the molecule (top panel) with ei = 0, £2 = —0.329, 
V = 0.0109, 7n,H = 0.03, all in units of eV, and a = 0.55^, 
T = 298 K. 7 d = 0, 0.01, 0.1, 0.3 eV. 


III. RESULTS 

A. Two-state diodes 

The functionality of the junction as a diode is demon¬ 
strated in Fig. [2] Inelastic-environmental effects deteri¬ 
orate the rectifying operation, as expected, yet the junc¬ 
tion rectifies substantially even when 7 ^ is order of the 
tunneling energy v. When inelastic effects are strong, 
7 d ^ V, the junction supports only i? ~ 50. Naturally, 
these values reduce when £2 is placed closer to the Fermi 
energy, when the gap ei — £2 is reduced, and when 'yL,R 
takes larger values, as was demonstrated in Ref—. 

To understand the decline in the diode operation un¬ 
der environmental effects we display in Fig. [3] the energy 
resolved currents at the two terminals, /^(e) and Inie). 
These currents are calculated from the integrand in Eq. 
([ 6 ]) and the corresponding equation for the right terminal. 
We find that at small biases (0.05 V) the resolved current 
arranges a single peak located close to the Fermi energy. 
The role of the environment is to broaden this peak, and 
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FIG. 3. Energy resolved currents in the coherent limit = 0 (black) and under inelastic effects 7 d = 0.05 eV. In the latter 
case \lL{t)\ 7 ^ |7_R(e)|, i.e. the current profile of electrons eirtering the junction is different from the one leaving it. The three 
panels follow the configurations of Fig. [T] (a) Low bias, when inelastic effects play a little role, (b) At peak voltage ~ 0.6 V the 
two states are in resonance, ei + aeE/2 ~ £2 — aeF/2 = —0.165 eV. (c) At large negative voltages —0.6 V, ei + aeF/2=0.165 
eV and £2 — aeE/2 = —0.49 eV. In panels (b)-(c) we further depict the resolved current under the dephasing probe condition 
when energy relaxation is not allowed (dashed line). Parameters are the same as in Fig. [2] 


to slightly enhance the current due to the contribution 
of inelastic terms. The behavior of the resolved currents 
at the voltage 0.6 V is displayed in Fig. [3Kb). This volt¬ 
age arranges a resonance conhguration most benehcial 
for conductance, see Fig. HKb). When = 0, a single- 
large contribution to the current shows up (clipped here 
for clarity), located around e = -0.165 eV, corresponding 
to the energy of the two molecular sites. When we allow 
energy exchange with the environment, electrons with 
energies reaching = 0.3 eV contribute, see |/_R(e)|. 
These electrons dissipate their energy on the junction, 
and are re-emitted at the other terminal with lower ener¬ 
gies, down to = —0.3 eV, see the behavior of |/L(e)|- 
Under the reversed bias —0.6 V, corresponding to the 
conhguration of Fig. He), an analogous behavior takes 
place, see Fig. He). 

We further present in Figs. |SKb)-(c) the resolved cur¬ 
rent under the dephasing probe condition, when decoher¬ 
ence effects take place without relaxation. In this case, 
energy conservation is enforced and electrons on the junc¬ 
tion do not absorb/dissipate energy from/to the thermal 
bath, thus |/L(e)| = |/_R(e)|. The effect of the environ¬ 
ment is to broaden the coherent current prohle, again dis¬ 
turbing the operation of the diode. Interestingly, the mi¬ 
croscopic differences between the two probes as reflected 
in Fig. Odo not translate to the macroscopic operation of 
the device in the present model, and the total-integrated 
currents under either probes are almost identical. This 
observation however is not general, and deviations be¬ 
tween probes (emulating different microscopic processes) 
may clearly affect functionality, particularly in multi-site 
chains. 

The departure from the strict tunneling picture limits 
the operation of the junction as a diode since the contrast 
between resonant transmission and tunneling no longer 
stringently holds with inelastic electrons contributing to 


the current. Note that 71 ,,ij, the coupling of the molecule 
to the metals, corresponds to levels’ broadening, and it 
affects the diode in a parallel manner as the tunneling 
current becomes over-dominated by ballistic electrons. 
This effect is presented in Fig. [H displaying contour plots 
(log scale) of rectification ratios as a function of '^l,r and 
applied voltage bias at different values for ■jd- 

In Fig. |S]we further present the rectification ratio at 
peak voltage, R{Vp), as a function of 7 ^. Weakly hy¬ 
bridized junctions operate better as tunneling diodes, yet 
they are more susceptible to inelastic effects. When the 
levels are largely broadened with 'yL,R = 0.1 eV, the junc¬ 
tion still rectifies well with R ^ 100 in the coherent limit, 
and it is resistant to the environment as long as 7 d O.I 
eV. 


B. Four-state diodes: sites with quasidegenerate states 

The rectifying behavior of a tunneling diode could be 
enhanced signihcantly if we mange to attenuate the tun¬ 
neling current under reverse bias, while keeping intact the 
resonant current in the forward direction. We recall that 
in an A^-site molecular chain the tunneling conductance 
decays exponentially with distance, while the resonant- 
coherent current is independent of chain length^ (leaving 
aside decoherence and dissipation effects). We thus sug¬ 
gest to generalize the two-state tunneling diode of Fig. [T] 
and include four molecular states, two sites within each 
conjugated unit, see Fig. [51 In the energy representation, 
each conjugated unit includes two closely spaced orbitals 
(at zero bias) instead of a single level. This principle was 
employed in RefJ^ to experimentally realize highly rec¬ 
tifying junctions, albeit employing only one conjugated 
unit and a long saturated tail. The Hamiltonian of the 
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with Hi 2 describing the conjugated units and V 12 and 
V 21 comprising the tunneling terms through the cr bridge, 


= 


€i + aeF/2 vi 

vi ei + aeVl2 


; V 12 = 


0 0 
V 0 


fe2 - aeVl2 V 2 ^ t > ^ 

^ 2 =' V2 e2-aeF/2]’ = ( ^ g 


The two orbitals within each conjugated unit are made 
close in energy; ui ^2 determines departure from degener¬ 
acy. 

To make a meaningful comparison to the two-site junc¬ 
tion, we employ here the same parameters as in Fig. [2] 
The IV characteristic of the model is displayed in Fig. 
ini demonstrating a very large rectification ratio in the 
coherent limit, but a strong decline in operation under 
environmental effects. The four-site design thus offers 
an improved diode functionality, yet a greater sensitivity 
to thermal effects since a minuscule tunneling current is 
more easily outplayed by inelastic processes. 


FIG. 4. Map of rectification ratios as a function of applied 
bias and metal-molecule hybridization, (a) Coherent limit, 
(b) 7 d = 10 meV, and (c) 7 d = 100 meV. Parameters are the 
same as in Fig. [2] 



FIG. 5. Decline of rectifying behavior due to environmental 
interactions with the parameters of Fig. [2] Filled symbols at 
the left boundary mark corresponding rectification ratios in 
the coherent 7 ^ = 0 limit. 


four-site molecule takes the form 



IV. CONCLUSIONS 

We analyzed the operation of molecular tunneling 
diodes comprising conjugated-saturated-conjugated seg¬ 
ments under environmental effects inducing decoher¬ 
ence and energy relaxation. Specifically, we employed 
a two-site molecular junction with parameters corre¬ 
sponding to the I,2-bis(4-(phenylethynyl)phenyl)ethane 
junctional. This system operates as a diode since (i) an 
asymmetry is built into the molecular backbone, and (ii) 
an internal potential drop develops within the molecule 
in response to applied bias. For these reasons, we ex¬ 
pect our findings are fully general to constructions with 
different backbones or system energetics, provided they 
follow these design principles. The contrast between res¬ 
onant and tunneling conductances is the key to the op¬ 
timization of the diode, with rectification ratios reaching 
three orders of magnitude in coherent situations. The 
diode performs well under environmental effects as long 
as 7 d < V, providing R ~ 100 — 500. The decline in the 
operation of the diode under thermal effects ensues from 
the contribution of incoherent electrons, playing down 
the contrast between resonant transmission and tunnel¬ 
ing conductance. We further proposed the design of an 
improved rectifier, suggesting a chemical backbone sup¬ 
porting two closely spaced orbitals within each conju¬ 
gated unit. This structure provides R ^ 10® within our 
parameters when 7 d = 0, though this type of diode shows 
a greater susceptibility to environmental effects. 

An inelastic vibration-assisted diode obeying a similar 
principle of level alignment was recently investigated in 
Refs4^i^, with a weak, phonon-assisted hopping energy 
instead of the direct tunneling v. Despite the different 
conduction mechanism (resonant vs. hopping), this setup 
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Bias (V) 


FIG. 6. Top: Scheme of a four-site diode in the site repre¬ 
sentation (left), and energy representation within each con¬ 
jugated moiety (right). Bottom: Operation of the junction 
as a diode (a) IV characteristics, (b) rectification ratio, and 
(c) level diagram with (dashed lines) chemical potentials and 
(full lines) the four eigenenergies of the molecule. Parameters 
are ei = 0, £2 = —0.329, v = 0.0109, 7 n,H = 0.03 eV, and 
a = 0.55 as in Fig. [2] as well as vi = V 2 = 0.03 eV for the 
tunneling energy within each conjugated unit. The temper¬ 
ature is set a,t T — 298 K and 7d was varied as indicated in 
the figure. 


produced IV characteristics closely resembling the curve 
in Fig. EI 

The probe technique employed here for mimicking de¬ 
coherence and inelastic effects can be feasibly imple¬ 
mented within DFT -h NEGF calculations to receive 
first indications on the role of thermal effects on electron 
transport in nanojunctions^i^. Our work here is sig¬ 
nificant beyond the specific analysis of tunneling diodes 
as we demonstrate for the first time the implementa¬ 
tion and utility of the voltage probe method to molec¬ 
ular electronic problems under large biases, introducing 
inelastic effects into far-from-equilibrium molecular con¬ 
duction. Our approach allows the interrogation of the 
role of the environment on other elemental molecular 
electronic functionalities including negative differential 
conductance^, single-molecule transistors^, and nonlin¬ 
ear, phonon-assisted, thermoelectric energy conversion^. 
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